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nils  Hemorandun  describes  one  of  tiie  first  Project  RAND  studies 
in  a  continuinc  research  program  on  heat  and  mass  transfer  in  boundary- 
layer  flows.  Qiis  vrork  formed  the  basis  for  subsequent,  more  cosqplete 
studies  on  the  problem  of  re-entry  heating,  which  were  summarized  in 
RM-2516,  A  Review  of  Bi«*'»*y  Boundary  Layer  Characteristics,  by  J.  P. 
Gross,  J.  P.  Hartnett,  D.  J.  Masson,  and  C.  Gazley,  Jr.  Althou^  the 
work  was  done  in  1956,  major  portions  of  the  analysis  still  represent 
a  unique  contribution. 

!Qxe  results  reported  herein  were  first  presented  (orally  by  C. 
Gazley,  Jr.)  at  the  Mass-Xransfer  Cooling  Synposium  held  at  Oie  RAND 
Corporation  in  June  1957* 


SUMMARY 


This  Memoraiidum  reviews  the  jroblen  of  a  lamloar  boundary  layer 
on  a  flat  plate  with  mass  transfer.  Previous  work  Is  reviewed  cuad 
the  fundamental  equations  are  derived.  The  ease  of  an  Incompressible 
laminar  boundary  layer  without  heat  transfer  with  hydrogen,  carbon 
dioxide,  and  Iodine  Injection  Is  solved.  Velocity  and  concentration 
profiles  and  skln-frlctlon  coefficients  are  calculated  and  discussed. 
It  Is  shown  that  surface  Injection  of  a  foreign  material  Into  the 
boundary  layer  reduces  the  skln-frlctlon  coefficient.  The  stability 
of  an  Incompressible  laminar  boundary- layer  with  mass  transfer  Is 
estimated  using  the  Lin  approximation  and  is  shown  to  decrease  as 
the  molecular  weight  of  the  Injected  siibstance  decreases. 
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I.  unsoDiXTnm 


Ihe  binary  lanlnitf  boundary  layer  may  be  defined  an  a  boundary 

layer  In  vhlch  a  foreign  subatanee  la  dlffualng.  For  many  yeara  aelentlata 
connected  vlth  the  aircraft  Induatry  hanre  been  aeardilng  for  poaalbla  methoda 
of  controlling  the  boundary  layer  on  the  aurfacea  of  aircraft  and  mlaallea. 

It  haa  been  shoun  experiaantally  that  the  irietloin  coefficient  for  the  turbu¬ 
lent  boutidaxy  layer  nay  be  trcn  two  to  five  tlmaa  aa  great  aa  that  for  a 
laminar  layer.  IherefOre,  aueeeaaful  efforta  to  control  the  transition 
point  from  to  turbulent  boundary-layer  flow  would  be  rewarded  by 

elgtilflcant  aavlnga  In  power  requlraaMnta  for  the  vehicle,  ibia  has  re¬ 
sulted  In  SOBM  critical  analyses  of  the  stability  of  the  boundary 

layer,  l.e.,  the  causes  for  the  Instability  and  the  conditions  under  vhlOh 
It  occurs. 

Several  methods  for  controlling  the  boundary  layer  have  been  proposed. 
Raaiovlxig  pert  of  the  boundary  layer  throu^  a  slotted  or  porous  surface 
prevents  thickening  of  the  boundary  layer  and  consequent  turbulence.  How¬ 
ever,  this  iBBthod  also  increases  the  friction  coefficient  and  thus  Is 
limited  In  Its  applicability.  It  Is  also  possible  to  design  the  airfoil 
In  such  a  uay  that  transition  Is  delayed  as  long  as  possible.  However, 
structural,  axnor,  and  fuel  requirements  play  a  very  inportant  part  In  the 
design  of  an  airfoil,  axid  It  may  be  necessary  to  abandon  some  conditions 
In  order  to  fulfill  others.  Beyond  that,  the  difficulties  In  aerodynamic 
design  at  svqpersonlc  speeds  are  wall  known.  Finally,  It  has  been  shown 
that  the  boundary  layer  may  be  stabilized  by  cooling.  This  presents  three 
possibilities:  (1)  direct  removal  of  heat  by  conductions  to  a  heat  sink, 

(2)  Injection  Into  the  boundary  layer  of  a  Ughtwel^t,  hl^-heat-capaclty 


gas,  or  (3)  use  of  a  solid  surface  vhiOh  would  sublime  or  vaporize  into  the 
bouxtdary  layer  at  elevated  tergQperatures.  The  latter  two  possibilities  nay 
Involve  a  binary  Inmlxvar  boundary  layer. 

Gbe  direct  rsnorval  of  heat  by  transfer  throuf^  the  skin  to  a  heat  sink 
can  be  shown  to  be  an  order  of  oa^iltude  less  efficient  than  the  latter 
two  cooling  methods,  m  the  ease  of  extrsnely  hl^  fllf^t  speeds,  a  systen 
enqitloytng  this  method  ml^t  become  incapable  of  handling  the  heat  loads  al¬ 
together.  Eckert,  Schneider,  and  K&hler^^^  have  shown  the  beneficial  re¬ 
sults  obtained  by  Injecting  a  U^tweight,  hlgh-heat-eapaclty  gas  into  the 
boundary  layer.  Obey  do  not  take  into  account  the  fact  that  hl|^  rates  of 
injection  of  a  llc^tweii^t  gas  act  in  a  destabilizing  manner  and  mi^t 
bring  about  early  transition  to  turbulent  flow.  As  pointed  out  above,  this 
is  undesirable  from  a  design  point  of  view.  Injectlxig  materials  of  hl#er 
molecular  weleiit  into  the  boundary  layer  will  not  give  the  low  friction 
coefficients  that  are  obtained  with  hydrogen.  However,  the  delay  In  tran¬ 
sition  to  turbulent  flow  may  well  be  the  deciding  factor  here.  A  quali¬ 
tative  stability  analysis  will  be  neoessaxy  to  Indloate  the  paroper  cooling 
method. 

A  crude  method  for  detemlnlng  the  stability  of  a  compressible  laminar 
flow  is  to  observe  the  varlatlen  of  the  quantity  pdu/dy  throu^  the  bounda¬ 
ry  layer.  If  this  quantity  has  a  point  of  infleetian  for  some  u/u  >  1  -  1/M 
(see  Curve  1  of  text  figure  below),  the  flow  is  unstable  at  eufflclently 
hl^  Reynolds  nuabers.'  *  Curve  2  Bhowe  the  eases  for  a  stable  flow  and 
Curve  3  Indicates  a  neutral  sltuatlen.  If,  then,  the  shear  remains  essenti¬ 
ally  constant  for  the  addition  of  small  amounts  of  forelga  material  to  the 
boundary  layer.  It  might  be  possible  to  Improve  the  shape  of  the  pdu/dy 
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ouxve  by  Increasing  the  density  and  decreasing  the  viscosity  of  the  binazy 
boundary  layer.  It  should  be  noted  that  this  alone  vlU  not  bring  about 
the  desired  result,  for  It  Is  essentially  a  change  of  sign,  and  not  magnitude 
only,  that  Is  of  Interest  h«:«.  Previous  vork  has  shown  that  the  Intro¬ 
duction  of  a  forelfp,  diffusing  substance  Into  the  boundary  layer  produoes 
a  change  In  the  velocity  profile  by  the  addition  of  a  vertical  velocity 
component.  Cie  magilfade  of  this  effect  will  depend  on  the  porapertles  of 
the  Injected  material.  In  general,  low-density  substances  will  be  more 
destabilizing  because  they  are  more  effective  In  reducing  the  flow  acceler¬ 
ation  near  the  wall.  !Bils  Is^pUes  the  appearanoe  of  an  Inflection  point. 

Hlj^ier-denslty  matwlals  ^ch  do  not  decelerate  the  flew  significantly 
will  as  a  consequence  reduce  the  skin  friction  by  only  a  small  amount, 
niese  two  bpposlte  effects  must  than  be  considered.  A  Im^^denslty  material 
will  cause  a  large  reduction  in  skin  ITictlen  but  may  destabilise  the  boundary 


layar  to  such  an  extent  that  transition  vUl  occur*  A  hl^-denslty  Injection 
naterlal  vlll  not  bring  about  an  early  trensltlonf  but  the  resulting  akin* 
inetlon  reduction  vlll  be  minimal  *  Sane  congircnlse  vUl  be  necessary. 

On  the  basis  of  the  ahcnre  arguments^  It  vould  seem  reasonable  to  search 
for  aaterlals  possessing  relatively  high  densities  and  low  viscosities*  Tuo 
such  substances  'lAil^  are  feMlble  iron  a  practical  standpoint  are  carbon 
dioxide  and  Iodine.  Both  of  these  aaterlals  sublime,  so  that  it  may  be 
possible  to  employ  them  vlthout  the  attendant  use  of  trexisplratlon-cooUng 
plumbing.  A  qualitative  analysis  vlll  be  a^plJQyed  using  the  Isotbemal 
boundary  layer* 
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n«  equations 

Qie  equations  governing  a  laminar  binary  boundary  layer  flav  vlU  be 
derived  borlefly  foUovlng  oloaely  the  analysis  of  Ball*''''  la  general,  ve 
shall  be  eoneexned  vlth  a  aeiltloonpooent,  nonreaetlng  aystenu  Our  furtiier 
assuqptlons  vlU  be 

!•  Isotropic  aedltm 

2.  Single-phase  and  locally  homogeneous  aedluB 

3.  Ho  external  force  field 

Transport  phencmena  linearly  dependent  iqpan  property  gradients 
3*  Diffusion  velocities  small  In  absolute  magnitude 
6.  Stea^  state 

Since  the  system  Is  assumed  to  be  xionreactlng,  the  generation  term 
does  not  appear  and  the  continuity  of  mass  for  each  coeiponent  may  be  ex¬ 
pressed  as 

^  (1) 

vbere  n^  Is  the  number  of  molecules  of  species  1  per  unit  volume  and 
Is  the  diffusion  velocity.  Sunalng  all  the  coBQ>onents  gives 

^  +  7  •  =  0 

vhere 

» -  Z  “i“i 

1  1 

A  ,A  A  , 

PV  m^n^(V  +  Vj) 

1 


(2) 

(3) 

W 


and 
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the  last  equation  defining  the  oenteavof-oass  velocity  V.  ^e  diffusion 
mass  velocity  is  measured  hy  the  relative  mass  motion  of  a  pertiexilar 
ocnponent  vith  respect  to  the  center-of-mass  velocity: 


ii  «  n^nj^^  »  m^n^(7^  -  ^ 


(5) 


fOierefore 


A 

J4 


A  ^ 
V  + 


(6) 


and^  finally 


dc, 

dt 


+  V  •  =  0 


(7) 


Equation  (7)  is  referred  to  hereafter  as  the  basic  diffusion  equation. 

Ihe  equations  of  motion  are  simply  the  veil-known  Haviar-Stokes 
equations.  Which  have  been  derived  in  detail  in  Ref.  4: 


p  ^  -  -Vp  +  V  •  (K  -  I  u)  •  ^  +  u  (7  V  ♦  V*  ^ 


(8) 


It  should  be  noted  that  all  property  values  are  taken  as  the  looal  values 
of  the  mixtures. 

The  energy  equation  can  be  obtained  by  considering  the  time  rate  of 
change  of  energy  content,  the  convection  of  energy,  and  the  flow  of  work. 
Kegleetlng  dlfHerenoes  between  co^onent  kinetic  energies  gives 


ihere 


i 


viscous  stress  tensor 
energy  flow  vector 


(9) 
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The  second  term  In  tiie  braclcets  describes  the  diffusion  of  work  due  to  the 

action  of  the  stress  field;  the  last  term  describes  the  diffusion  of  energy 

rather  than  heat.  In  a  multicomponent  system  careful  distinction  should  be 

A  ^ 

made  between  heat  flow  q  and  energy  flow  Q,  vdilch  includes  the  effect  of  the 
diffusing  components. 

Introducing  the  equations  of  motion  and  continuity  results  in  the 
following  simplification: 


A  A  A  Z  A 

pV  •  Vh  “  V  •  Vp  +  V  •  q  -  T  »  TV  ■  0 


(10) 


Now,  it  may  be  shown  in  classical  thermodynamics  that  the  rate  of 
change  of  enthalpy  is  equal  to  the  change  of  heat  occurring  throu^  any 
thermal  transfer  mechanism  (closed  system)  plus  the  heat  change  brought 
about  by  concentration  variations  (open  system) 


„  f!i  .i  Je 

dt  *  dt  V  ^1  dt  p  dt 


(n) 


idiere 


.th 


*^1  >  chemical  potential  of  1  ccnponent 
S  «  entropy  of  the  system 

In  order  to  solve  Eq.  (ll)«  the  first  two  terms  on  the  rl^t  must  be 
evaluated  or  e3q)ressed  more  directly  in  the  properties  of  the  system.  In 
irreversible  tbexuodynamics^^^  it  has  been  shown  that  an  entropy  continuity 
equation  may  be  written: 


where 


da  ..  ^ 

P  ^  «  •  <7 


s  ■  entropy  flux  vector 


(12) 
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a  ■  rate  of  Irreversible  entropy  production 
The  texn  on  the  left  Indicates  the  total  change  of  entropy  within  the 
system.  V  •  s  Is  the  reversible  flow  of  entaropy  thzou^  the  bounding  B\ir> 
face  and  a  1b  the  rate  of  irreversible  entropy  production. 

The  reversible  entropy  flow  vector  Is  given  by 

Ts  -  Q  -  Z  4*  J.  (13) 

1  ^  ^ 

and 

Q  •  q  +  Zi  h.  J.  (14) 

1  ^ 

\diere 

>  partial  enthalpy  of  1^  component 
Combining  Eqs.  (lO),  (ll)>  (12) >  and  (13)  we  obtain  the  following  ex* 
pression  for  the  rate  of  Irreversible  entropy  production: 

Tb--q  .  VlnT-Ij.  .X.  ♦T.w  (15) 

1  ^  ^ 


where 


\  -  TV  (-i)  +  ^  ▼  In  T 


(16) 


(7^ 

The  lElnetlc  theory  of  dilute  gases' ' '  gives  eaqpresslons  for  the  Hux 
vectors  for  energy  and  diffusion: 


-  X*  VT  -  p  I 


J  ^ 


a 

q 


(17) 
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"  ■  4“  ^  r  I 

i 

T 

The  are  gsoerallted  dlfAiaioo  coefficients^  the  are  generalized 
tfaemal  diffusion  coefficients^  and  X*  Is  Ibe  ’^st"  liieziial  conduetlvlty** 
HSie  Is  given  as 

^  M  fe>)  ’"J  *  <Vl  -  '!>’  ^  *  <19> 

j  \  J/t,p 


vhere 

>  the  partial  molal  volume 
lOr  a  mixture  of  Ideal  gases^  Bq.  (19)  becomes 

^1  “  r  ^  ^  P  (19a) 

Fot  a  perfect  gas,  Eq«  (19a)  reduces  to 

-  R^T  ▼  In  p  +  R^T  7  In  (-|)  (I9b) 

•X  ,  »  V  -X  X.  A 

Substituting  d^  Into  Eq.  (19)  and  noting  /L  d^  >  0  results  In  d^  >  X^» 

Equations  (17)  and  (iB)  Indicate  a  linear  relationship  between  the  fluxes 

q  and  ^  and  the  "Anrces"  7  In  T  and  !Die  Onsager  reciprocity  relations 

which  are  valid  for  such  linear  relations  give 

"I  -  ’>]  ®U  ■  ®J1 


*  !Ihl8  Is  the  thezmal  conductivity  whldi  describes  heat  transfer  sole¬ 
ly  by  conduction  without  the  convective  effects  of  the  Induced  thenaodlffUslon 

flov. 
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Finally^  using  Eq. 


(4)  and  el  littlnatlng  ftoa  Eqs.  (I7)  and  (iS)  per¬ 


mits  us  to  write 


4 


(21) 


and 


^  -  d^)  -  V  In  t]  (22) 


T 

idiere  k  Is  the  usual  thermal  conductivity  and  k  Is  the  thermal  diffusion 
coefficient.  It  should  be  noted  that  k  includes  the  conducting  effects  of 
the  induced  diffusion  stresn  and  does  not  reduce  to  X  even  for  ^  ■  0. 


q  .  -x'vT'S  -  m  (22e) 

t 

It  mie^t  be  noted  that  k  and  X  differ  by  a  negligibly  ■nail  aeount  for 
materials  with  hl#i  molecular  weieJits.  If  we  are  concerned  with  li(^t  perti- 


cles  such  as  electrons,  then  this  difference  becomes  significant. 

Since  the  problem  concerns  a  two-cosqponent  system,  the  eolations  will 
be  simplified  for  this  epeclal  case.  Oie  components  are  related  by  the 


following  cenditions: 

®1 

^1  *  ^2  "2 

Oien  the  eqoations  can  be  put  in  the  fOxm 

It  ® 

P  ®  ^  -  0 


(23) 


(2) 


(24) 
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vhere 


P^t  “  “''^P  +  V.(k-|  |i)t7  .  V  +  u(VV  +  (8) 

«  ■  •  <*®1  (25) 


-,n 


'12 


7e^  ♦ 


Mg-M, 

“'  k 


(1 


Cj^)V  In  p  +  ac^il 


a 
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m.  BOUNDAHY-IAYBR  BCU4TI0HS 


Consider  the  problem  of  a  binary  boundary  layer  on  a  flat  plate.  Die 
equations  gonreming  the  flov  of  a  single  cooqponent  gas  over  a  flat  plate 
are  veil  known  and  have  been  extensively  investigated.  A  solution  of  the 
binary-boundary-layer  equations  will  require  a  careful  order-of-magnltude 
analysis  to  determine  if  any  sijapllflcatlons  are  permissible.  Ibis  has 
been  done  by  Hall.^  ^  The  analysis  of  the  mass-diffusion  will  be  repeated 
here  because  of  its  importance  in  the  problem. 

Let  two  new  variables  X  and  Y  be  defined  axid  nondimenslonallze  u,  v, 
and  c  as  follows: 

L  ■  characteristic  dimension  (x-direction) 

5  a  thickness  of  boundary  layer 
u  ■  veloci-ty  in  the  x-dlrectlon 
V  a  velocity  in  the  y-directlon 
c^^  a  concentration  at  the  wall 
e  a  properties  referred  to  edge  of  boundary  layer 
Substitution  of  new  variables  into  Eq.  (2^)  results  in  the  transfor¬ 
mation  r  _  f  _  ^ 

“e^  Ir] 

n  1  2-  J  -  *  ”2  ’  7/1  J-P  0  *  .  a  In  T 

°i2  p  ax  I  p  [ax  Mp^  ®i'^  "  ®r®iw  ax  i'^  "  r  ax 

(27) 
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Hhe  two  tezms  on  the  left  side  of  Bq.  (22)  will  be  coBqpered  first*  An 
order-of-nagnltude  axuJysls  of  the  continuity  equation  Indicated  that 

\X 

V  >  0(v  u).  Dlls  hnpUes  that  both  terns  are  of  0(-r).  Since  It  Is 

postulated  that  L  »  b,  only  the  first  tern  on  the  rlj^t  side  of  0(-^)  will 

be  Isqtortant  under  the  conditions  stated.  If  a  diffusion  boundary  layer  Is 

to  exlstj  both  sides  must  be  of  similar  order  of  magnitude: 
u 

0(j^)  ■ 

This  may  be  rewritten  as 

o(^)  -  0(4-)  (89) 

L  b'^Sc 


If  the  diffusion  boundary  layer  Is  to  be  similar  to  the  momentum  boundary 
layer^  then  It  Is  evident  that 

sc  -  0(1)  (30) 

nothing  further  can  be  stated  concerning  the  terms  within  the  brackets 
until  more  Information  Is  available  about  the  derivatives  and  their  coef¬ 
ficients.  Consequently^  the  final  form  of  the  dlffUslon-boundary-layer 
equation  Is 


pu 


“if- 


(31) 


^e  other  boundary-layer  eijuatlons  are 


ax  7sy 


(32) 


(33) 
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2l£ 

ay 


«vli 


'V 


^  - 

ay 


a 

?y 


(3lv) 

(35) 


-  »  ^  ”*“12^  ’  l>®12<%  ■  V  »  ■* 

h  ■  -  'l'  |y  ^  ^  <36) 

If  WB  assume  that  the  similarity  coodltlon  applies,  l.e.,  that  the  depeod- 
enee  upon  x  and  y  can  be  eonbined  Into  a  single  variable  giving  the  normal 
position  vlth  a  scale  adjustment  according  to  the  position  along  the  plate, 
then  the  system  of  partial  differential  equations  may  be  transformed  Into 
a  set  of  oordlnary  differential  eciiatlons.  The  variable  T)  Is  defined  as  a 
dimensionless  (fiantlty  by 


T1 


(37) 


It  should  be  noted  that  the  Introduction  of  this  parameter  may  restrain  the 
functional  dependence  of  the  unknowns.  For  instance.  It  will  be  shown 
later  than  v(x)  must  be  proportional  to  1/  -/x^ If  tbo  equations  are  to  exhi¬ 
bit  unidependency.  For  the  case  of  sublimation  or  evaporation,  the  species- 
conservation  eqioatlon  at  the  sxurfhce  yields  the  required  v(z)  proportional¬ 
ity.  Firthemore  we  define  a  stream  function 

t  (38) 

ihldi  satisfies  the  continuity  equation 

Pefy  •  Pefx  "  ’P  <39) 
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Zf  -these  ftine-tlons  ere  aubs-ti-tu-ted  into  Ihe  boundary-la^  eguatlons,  -they 
take  the  A»m 


*lfi  (f' [S;  Is  ■  ®  * 

" Is  ^ ^ * ‘‘pj ■ 


Becapuse  -the  system  of  equations  Is  of  the  seventh  order.  It  vUl  be  neoes- 
saxy  to  provide  seven  boundary  conditions.  Siese  will  vary  vith  the  -type 
of  mass  transfer  that  ve  are  interested  in.  Eckert^  Schneider,  and  Kfihler 
have  solved  the  problaai  for  the  ease  of  injection  of  hydrogen  into  a  stream 
of  air.^^^  At  -the  outer  limits  of  the  boundary  layer,  the  condi-tionB  vUl 


alveys  be  -the 


y  -  e 


Ih  terms  of  our  nev  variables,  this  appears  as 


71  -  « 


T  .  T 


Ohese  eolations  have  been  solved  by  Eckert  for  -the  ease  of  a  porous  plate. 

He  found  -that  by  injecting  a  aautll  aanunt  of  li^tueight,  hi^-heat-eapamity 
gas,  it  is  possible  to  reduce  considerably  -the  aamunt  of  heat  removal  from 
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Um  ]^tt.  Bis  sfastlona  wsrs  sjtopliflsd  In  tbs  ssnss  that  hs  assuasd  tbs 
tbsznal  diffusion  oosfflelsnt  to  bs  ssm.  Hs  did  not  ixnrsstigsts  tbs  sts^ 
bUlty  diarsetsristies  of  tbs  systm, 

Sis  ooBg>lsidty  of  tbs  system  of  e<|iatlQg&s  ubidi  must  bo  solved  In  order 
to  ijonrestlsats  tbs  velocity^  tanperaturs,  and  oonoentration  profiles  is  such 
that  it  is  neesssaxy  to  perform  a  madiine  calculation*  Firsts  tbs  sipiatlans 
must  be  put  into  a  systaa  of  first-order  Runga-Kutta  e<|Aatioos.  We  let 

f  -  Yq  ®  "  ^3  ®1  “  ^5 

t 

Z"  ■  yi  t’  -  y,^  ®1  ■  ^6  ^^5) 

p 

t  III  III 

-  yg  T  -  y^  ®i  •  ^6 

N  ■  i-; 

Sis  boundary-layer  squatians  may  be  put  in  tbe  fOUoving  forms: 


*  ®1 

Uhere 

•  1  (46b) 

Ai2  -  0  (46o) 

-  0  (46d) 

A21  -  0  (46s) 

^22  ■  “  yj) 

3 

-  1  (46g) 
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*31  -  ® 


*  ^  «  V  vl 

*33  ■  -1^ 

\  -  -  — 

|i 


-  -V6 


r  yji  •  .  X  ^4  ^4 


y?l 

-  oQTjCi  -  yj)  “5  > 


H 


=^°p  V4 

P 


p  » 


f  „  1*  5[3 

|.  yg  -  Ojy^d  -  y,)  '  I??  ^ 

[{fe'"}  * '"5'*  -  ='3> 

-  Iwsd  -  yj)  ^  *  a'^jd  -  V  ^  •  ’^f6 


y* 


(46h) 


yj 

-  oor^d  -  yj)  ^ 


(46J) 

(l»6k) 


(4a) 


*  {^6  *  "  V 
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Sm  systen  of  o^tlonB  oontalned  In  Eqib.  (U6a  -  n)  appears  in  Runae-Kutta 
neliiod  as 

y'o  "  ^iP  H  “  ^6^ 

y[  -  yg/tr  yj  -  yg 

4  “  ^6^  ^6  “  '3(^0  •••  ^6^ 

y*  -  y,^ 

vhere  f^^  ^2*  and  f^  are  Obtained  ty  solving  the  system  of  Eqjs*  -  ■)• 
Before  solving,  howver.  It  Is  necessary  to  speeliy  four  other  boundary 
eonditlona.  At  the  aurfiaoe  there  ere  three  obvious  conditions: 
y  ■  0  u  ■  0 

T  -  (W) 

®1  “  ‘aw 

To  account  for  the  other  boundary  eoodltlan,  ve  observe  that  the  oonveetlve 
velocity  of  the  forel^x  gss  passing  throu^  the  tenuous  dumnels  of  the 
plate  vlU  be  large  compared  with  the  mass  velocity  of  the  air.  Ih  other 
words,  Mt  say  that  the  mass  of  air  ihl<h  enters  the  plate  la  negligible. 
Matheamtleally,  the  diffusion  current  for  the  air  Is  given  by 

%  - 

According  to  the  aseuqptlcn  that  the  nass  velocity  of  the  air  oust  be  zero 

for  y^O  (50) 

Tbe  surface  diffusion  current  far  the  forelgo  gas  becomes 

as 

Jy  ■  (^  •  ®lw^P¥%  ■  "  jy^w 


(51) 
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ZraziBlbinlJig  Into  our  now  Yarlables  gives 


It  diould  te  noted  that  liie  situation  Is  quite  different  for  a  subliming 
solid*  In  Ihls  ease  there  Is  no  question  that  the  mass  flow  of  air  back 
throuj^  the  plate  Is  zero.  It  Is  important,  hovever,  that  the  flov  of 
foreign  gas  Is  eoupled  to  the  temperature  boundary  layer  by  the  heat  of 
sublimation.  This  can  be  related  most  sliqply  through  the  Clapeyron-Clauslus 
eolation 


If  ue  are  oonsiderlng  eubUmation,  p  is  the  vagpor  pressure  over  the  solid; 
dV  a  Vg  •  AS  slnoe  >  V^;  and  Is  the  molar  heat  of  sublimation. 
Hov  if  ve  assume  that  the  gss  Is  Ideal  and  that  heat  of  sublimation  Is 
a  funetlon  of  temperature  only,  ue  obtain 

InPy  •  5  j  (5*>) 

there  Is  the  temperature  at  thloh  the  pressure  Is  1  aha. 

Ihe  Introduotlon  of  a  subliming  solid  coapllcates  the  boundary  con¬ 
ditions,  because  the  vail  tSBqperature  and  concentration  are  no  longer  vari¬ 
able.  For  the  case  of  transpiration  pooling  under  a  given  set  of  external 
flov  conditions,  the  mass  Injection  rate  (a  variable  boundary  Input)  controls 
the  concentration  of  the  forel^i  gas  at  ^e  surlhee  and  this  In  turn  per¬ 
mits  the  dioloe  of  a  variable  surface  tesperature  or  beat  transfer.  Sub¬ 
limation  cooling  Imposes  another  condition,  nasMly,  a  relationship  between 
the  surface  taqperaturs  azd  concentration.  The  heat  of  sublimation  fixes 
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dqpoDdBnoe  aB  ahova  t&  Eq.»  (5U).  Fmrthenaortt,  -ttie  hsat  traztsfar  at 
tlie  vail  la  ralatad  to  lha  heat  of  ahbllaatlea  in  Him  following  way  (aaauBdng 
DO  radiation  or  oonduotion  loaaea): 

It  ahould  be  noted  that  the  takaa  Into  account  the  laat  blocking  action 
of  the  foreiffx  gaa  dlffUalng  throui^  it.  Ibe  injection  rate,  therefore,  ia 
not  arbitrary  aa  in  the  preiviouB  eaae  but  ia  a  function  of  heat  Hux  and 
heat  of  aubUaatlon.  Chooalng  a  material  eaaantlally  lixea  the  heat  of 
Bubllaatian  (except  for  a  amall  dependence  on  tanperature)  •  A  vail  tanpe]> 
ature  nov  detemineB  the  vail  oonoentratlon,  hence  the  Injection  rate  and 
finally  the  heat  tranafar.  Thia  heat-tranafer  rate  nuat  aatlaiy  Bq*  (5^)» 
Ihe  problem  of  handling  the  complex  boundary  eondltiona  can  be  moat 
easily  taken  care  of  by  setting  up  another  set  of  linear  equations  Which 
must  be  solved  slaultanaously  with  the  original  ayatem.  The  RuDge->Kutta 
method  requlrea  only  initial  conditions,  and  so  a  series  of  Initial  guesses 
as  to  the  original  conditions  will  be  required.  The  guesses  will  lead  to  a 
set  of  solutions  vhldi  may  then  be  compered  with  the  boundary  conditions  at 
the  outer  edge  of  the  boundary  layer.  A  perturbation  of  the  initial  guesses 
and  an  exeminatlon  of  the  results  of  the  perturbation  on  the  end  conditions 
vUl  then  permit  a  linear  coapensator  to  adjust  the  Initial  guesses  until  'ttie 
final  boundary  oonditloDs  are  met. 


iy»  BOUHDAICf-DtfHl  aEUBTiamf 


Ohe  solutloa  to  12ie  boundary-layor  aquations  hanrlng  bean  obtained.  It 
vlU  be  neoessaxy  to  Investigate  llie  stability  diaraeteristies  of  the 
Tni<njit»  ooBqnresslble  boundary  layer  vlth  subUnatlon.  Lin  and  Dunn  hsnre 
shoun  that  an  order-ofHBagnitude  analysis  of  the  oengplete  linearized  equa¬ 
tions  for  a  three-diaensional  disturbanoa  imposed  upon  a  tHD-dlmenslonal 
boundazy  layer  indicates  that  they  oan  be  reduoed  to  oudi  slaplf  texns  In 
the  first  approodaatlon.  Biese  simplified  equations  are  still  valid  at 
hl^  firee-strssB  Madi  nisid)ers.  Houever,  the  aeouraoy  of  a  first  approaci- 
nation  nay  not  be  great,  so  that  hif^ier-ardar  agproixinationo  nay  hsnre  to 
be  oonsidered.  lin  and  Dutm  also  found  that  the  eonelusion  readied  by 
Lin  and  Laes,  naiMly,  that  the  relationship  for  the  (haraoterlstic  values 
is  Independent  of  the  tea^erature  iluetuatlons  and  the  boundary  oonditions 
Isiposed  upon  than,  is  restrieted  in  ganeral  validity  to  subsonic  axid  slightly 
supersonic  Maoh  msabers*  At  hl^i  Mach  xnmbers  the  relation  for  the  oharac- 
terlstlc  values  depends  in  general  upon  the  thermal  boundazy  conditions. 
Furtheznore,  Lin  and  Dunn  have  developed  a  netttod  for  detezolnlng  the 
stability  charsieteristics  of  a  coaqiressible  boundary  layer  vhldi  -ttiey  fael 
is  quite  aeourate  mp  to  a  Madi  number  of  two.  It  is  believed  that  their 
results  can  be  extended  to  Macb-number  values  as  hi^  as  six. 

The  work  of  lin  and  Duzm  was  not  concerned  with  blnazy  compressible 
boundazy  layers.  Iheoretieally,  a  new  set  of  equations  must  be  obtained 
including  the  diffusion  equation.  The  mathsnatieal  difficulties  of  solving 
such  a  system  of  equations,  even  thou^  they  are  linearized,  is  very  great. 
Instead,  as  a  first  approodmatiom,  it  ml{^t  be  possible  to  use  the  velocity 
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and  temperature  profiles  obtained  from  the  solution  of  the  above 
equations  In  the  Lin  and  Dunn  method.  This  Involves  assuming  that 
the  concentration  at  the  siirface  Is  exceedingly  smedl,  that  the 
diffusion  velocity  Is  negligible,  and  that  the  effect  of  the  foreign 
material  In  the  boundary  l^r,  as  far  as  stabilization  Is  concerned, 
will  be  evident  only  on  the  basis  of  the  velocity  and  temperature 
profiles.  This  assumption  has  been  shown  to  be  correct  by  E.  E. 
Covert. 
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V.  HUMTOCAL  30UUTI0M 

Before  Eqs.  (U6a  •  m)  can  be  solved,  tt  will  be  necessary  to  Indicate 
some  relationship  between  the  thermal  properties  and  the  temperature  and 
concentration  at  any  given  point.  Ohe  boundary  conditions  must  be  more 
carefully  defined.  Eq^iation  (^2)  accurately  defines  the  relationship  at 
the  wall  between  the  velocity  (of  the  foreign  substance)  Issuing  from  the 
wall  and  the  concentration  of  the  foreign  substance  at  the  wall.  However, 
In  the  case  of  a  binary  boundary  layer,  the  temperature  difference  existing 
at  any  point  also  causes  a  diffusion  of  material.  Of  course,  this  Is  also 
true  at  the  surface,  and  therefore  the  boundary  condition  of  Eq.  ($2)  Is  no 
longer  quite  coorrect.  Instead,  it  can  be  shown  that  the  influence  of  the 
temperature  gradient  on  the  mass  transfer  at  the  wall  can  be  taken  Into 
account  by  Including  a  thermal*diffUslon  term  In  our  original  equation: 


Equation  ($4)  can  be  rearranged  so  that  It  Is  more  amenable  to  calculation. 
If  the  total  inressure  Is  p^  and  the  beat  of  sublimation  Is  a  function  (quad* 
ratio)  of  tesperature  only,  the  equation  can  be  put  in  the  fOocm 


In 


(Sj^  ♦  SgT  ♦  S 


>*1 


If  it  la  assumed  that  the  wall  concentration  will  always  remain 
Eq.  ($6)  can  be  integrated  and  rearranged  to  give 


(56) 


then 
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*“•!*■* 


(5T) 


Ihe  Mnm  'boundary  oondltiona  aay  aov  ba  taou^t  togattiar  wd  put  Into  tha 


nav  oocrdlnata  aystai  aa  IbUowat 

-  2 

yjC-)  -  0 

y3(-)  ■  1 

y^Co)  -  0 

7(0)  r  ^  y.  (0) 

-  S-,75JI  “  yjCff)  *  a(0)y5(0)  1 1  -  y5(0)j  ^ 


la  ^5(0)  ■  R 


ya(o) 


^  fy3(0) -y,(i)'| 

^  \  yj(^)y3ri)  j  yjiy 


+  -  y3(l)j 


in '‘•“2 


(58a) 

(58b) 

(58c) 

(58d) 

(58a) 

(58f) 


(58g) 


The  thermal  properties  must  nov  be  calculated  in  terms  of  the  temperature 
and  concentration  at  any  given  point.  For  ree^ans  indicated  elsevbere  in  this 
Memorandum,  the  diffusing  substances  chosen  vara  iodine  and  carbon  dioxide.  The 
thermal  properties  of  these  materials  were  investigated,  and  those  properties 
which  were  not  available  from  experimental  data  were  calculated  using  statis¬ 
tical  mechanical  formulas.  In  all  cues,  the  properties  were  represented 
in  a  power  form  rather  than  the  original  formula,  e.g.,  the  Sutherland  vis¬ 
cosity  formula.  In  this  way,  the  number  of  constants  required  was  reduced 
to  a  mlniiaum.  All  thermal  properties  have  been  normalized  with  respect  to 
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air  at  the  free-streaa  temperature.  Qieee  nonukllaed  properties  are  represented 
In  the  imionHlng  maimert 


Densltgr: 

0  ■  a  .  Viirj/yj 

(59) 

Viseosityt  k,  k. 

-  ^  ^5  .  V3  -  ^5) 

-  <h*  W 

(60) 

Schmidt  nuaher: 

_  ha 

*Be  * 

^3 

(61) 

■Biexmal  diffUelTltgr: 

“  ■  ‘W'3  * 

(62) 

Pure  oaeponent  ttiexmel  eepaeltys 

-  >^5  *  'W'3  *  ‘^T»3 

(63) 

Mixture  thenMl  espaoitgr: 

"p  ■  ‘ifl  *  V5  *  'W'a  *  *  ‘2^/5  *  V5»3 


Mixture  Hiexmal  oonduotlyitgr: 

* - - 

*1  ■  h  • 

1  •  y  y 

■  1  ♦  kjjj  — ^  Wg  ■  1  +  3ri^  ^^5) 

Mixture  noleoular  wight: 

"  Mj^  ♦  (Mg  - 


(66) 
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Sone  of  vbIuas  for  12x680  oonoteats  for  the  ease  of  Iodine  rapor  -  air  aai 
earbon  dioxide  -  air  miocturet  are  given  in  Table  1.  To  aolve  12ie  eqoatlona^ 
the  foUovlng  oonstanta  nuet  be  specified:  and 

fia  addition^  either  12xe  oonoentration  or  the  ta^pearature  at  12xe  vail  oust  be 
specified*  Ohe  gradlenta  of  'ttxe  tanperature  and  oonoentration  at  12ie  vail 
are  also  raqoired^  but  loaknoun*  It  vHl  be  necessary  to  guess  12ie8e  In  ae- 
oordanee  vlth  our  Iteration  system  so  that  the  defined  boundary  conditions 
at  12ie  outer  edge  of  12xe  boundary  layer  result  uQpon  Integration* 
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Table  1 

COlBTAiniS  FOR  MUISTlCOMFCaiail  FROFOtTY  BOPiffilORS 


Constant 

00^* 

WWg 

•air 

lodlne-alr 

! 

-0.341 

-0.8858 

C 

0.855 

0.615^ 

0.884® 

0.764 

1.15^ 

0.883® 

1.2405 

2.3456 

‘5* 

-0.2405 

-1.3458 

1.125 

1.125 

“t* 

0.625 

0.625 

1.2405 

1.0000 

-0.2405 

-0.67053 

^0 

1.022^ 

1.004® 

1.933^ 

1.613® 

1.76^' 

1.66® 

1.87'' 

1.7® 

hi 

0.63^ 

-l.is'’ 

0.11^ 

0.00® 

2.44^ 

-2.428^^ 

o.a® 

.0.768® 

•0.21 

-0.755 

hit 

0.071 

-0.071 

hr 

-0.0046 

0.0038 

0.877 

0.877 

^9 

-0.207 

-0.755 

*20 

0.121 

0.125 

’*21 

-0.0036 

-0.0048 

-0.075 

-0.075 

^23 

-0.0048 

0.0048 

1.6202 

0.5271 

^5 

-0.6563 

-0.0695 
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Constant 

OOg-alx 

Iodine-air 

1.15  X  10"® 

1.15  X  10“® 

0.859 

de 

6.07i»  X  10"3 

1.09  X  10“*'  1.56  X  10**' 

1.396 

1.1^  0.833® 

^30 

2.0976 

7.676 

1.5361 

0.481 

^alues  obtained  from  experimental  data. 
^Inensionlesa  abaolute  temperature  less  than  3. 
dimensionless  absolute  temperature  more  thAn  3. 
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VT.  ffWlBTTJWf  CgARAgEPiailCS 

Ab  indicated  porenrloualy^  12ie  staMUty  obaraoterlstlee  of  tbe  diffusion 
systMS  will  'be  exasdned  by  placing  tioe  tflqpearature  and  velodly  profiles 
obtained  firm  Hie  solution  of  Bq,*  (^7)  Into  the  Lin  axid  Duxm  neHiod. 

Die  Hieccy  proposed  by  Lin  and  Dunn  vlll  not  be  presented  here,  but 
raHiar  Hie  nmerical  procedure  to  be  used  in  calculating  Hie  neutral  sta^ 
blUty  curve,  first,  we  define  two  functions  v(c)  and  x(c)  as  fOUowa: 


Given  a  y^Co),  v(c)  and  x(e)  can  be  detemlned  for  any  y^^Ce).  Xhe  follow¬ 
ing  approxlnation  scheoia  is  suggested: 


.n+1 


;«1)  .  - liiAh - 

(69) 

(1  +  xu“)  ♦  X^^ 

f  (1  ♦  xu“)  +  xVl  XV^ 

(TO) 

[(1  ♦  X)(l  ♦  Xu“)J  l  +  xu“ 

In  Eqs.  (69)  and  (70),  0^(z)  axid  0^(z)  are  known  functions  of  z;  a  tabul¬ 
ated  set  of  Hiese  functions  may  be  found  in  Lin  and  Dunn's  report. Since 
X  is  usually  small,  initial  guesses  for  the  approodaations  sHime  may  'be 

0i(z)  -  ▼(«)  >1  “  0,(*)  (71) 
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for  a  given  value  of  7^(0),  v(e)  and  x(e)  may  Toe  detexmlned  for  a  given 
y^(e).  NOV  v(c)  defines  an  initial  z  \Aiich  la  used  to  obtain  an  initial 
value  of  u.  7rm  Eq.  (69),  la  calculated.  Hhe  and  are 

put  into  Eq.  (70)  to  calculate  a  new  m  '(Ms  vay^  it  is  possible  to 

obtain  a  solution  (z^u)  for  each  'value  of  y^(c).  Ghe  vave  number  is  nov 
calculated  from 


Ic 

o 


o  ^y^Co)  (i(o) 
^  ^^y75T 


B^yi(o)yg(o) 
‘  1,3(0)  „(o) 


+  k. 


11 


+ 


J 


(72a) 


vhere 


Vl  ■  [1  «  yi(c)/2]' 

[l-yi(c)/2j^ 


(T2b) 


and 

B-  (72®) 

and  k^  may  be  comquted  by  madiine  as  indicated  in  Ap¬ 
pendix  B.  Ibe  Reynolds  number  may  now  be  obtained  firon 


(73a) 


ilieore 


[yi(c)  -  y3^]y,(o)u 

^3 


1(c) 


(73b) 


vn«  SIMPUFIEP  SOLUTION 


Even  a  cursory  examination  of  the  equations  to  he  solved  vould  Indicate 
the  dlfflcxilty  Involved  In  making  reasonable  approximations  to  the  Initial 
values  of  the  system  which  the  Runge-Kutta  routine  requires.  Qie  stability 
of  the  Convergence  Scheme  was  not  studied;  but  It  may  be  presumed  that  If 
the  Initial  values  are  not  at  least  "approximately”  close  to  the  proper 
values,  then  It  nay  be  that  the  routine  is  nonconverging  and  will  never  give 
a  solution.  To  prevent  such  difficulties,  it  ml^t  be  Judicious  to  solve 
a  simplified  case  first  so  that  the  Initial  conditions  resulting  tim  this 
attesqpt  may  be  used  In  the  more  ccmplex  solution.  Time  limitations  also 
played  an  Important  pert  In  choosing  the  simplified  ease. 

If  we  presume  that  the  fluid  is  incompressible  and  that  no  heat  trans- 
fer  is  present,  the  equations  presented  prevloualy  reduce  to 


(7^) 

(75) 


!Ihe  boundary  conditions  for  this  fifth-order,  ordinary  system  of  differen¬ 


tial  equations  are 

f'  »  0 


^Iw 


1 


f  -  f. 


0  (wall  concentration) 
(injection  parameter) 


(76a) 

(76b) 

(76c) 
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To  UM  th«  Ruaifli-Katta  systan^  all  bou&dazy  eandltiotxa  vlU  h«ira  to  ba  apeel- 
flad  at  12m  vail,  l.a.,  yj^(O).  Thd  initial  condltlona  axv  deflnltaly  known, 
Itaa  valoeity  at  the  vail  la  zaro,  and  Iha  concantratlon  hare  may  ha  spaol* 
flad.  It  vUl  ba  necaaaaxy  to  fulfill  and  boundary  condltlona  such  Hiat 
concentration  and  the  concantratlon  gradient  bacone  zaro  at  the  outer  edge 
of  Hm  boundary  layer.  If,  using  Bqa.  (7^)  and  (75),  va  define 


^0  • 

f 

(77a) 

yi  • 

1/p 

(77b) 

^2  - 

J  d/dn  [l/F  df/dTlj 

(77c) 

®1 

(77d) 

tiien  the  ayatem  of  equations 

P  dCi/dl) 

appears  in  the  desired  linear  fom 

(77e) 

^0  • 

(78a) 

yi  - 

yj^ 

(78b) 

- 

(78c) 

yj  - 

yj/p 

(78d) 

y4  « 

vlth  the  boundary  conditions 

* 

(78e) 

yi(o) 

m  0 

(79a) 

y3(o) 

'  ®lw 

(79b) 

“  'v 

(79c) 

yi(e) 

»  2 

(79d) 

y3(e) 

»  0 

(79e) 

Knowing  y^iO)  and  y3(0),  find  it  necessary  to  choose  y^iO)  and  ygCO)  in 

such  a  way  that  the  last  tw  boundary  conditions  are  fulfilled. 

A  coirrer- 

gence  achame  was  used  vhich  depended  upon  a  grid  of  goeaaea  \iilcfa  vaa  than 
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latcrvolatsd  to  gLv*  the  proper  Initial  valusB.  Sia  solution  of  these  eqi- 
etiODs  (exsotly  slartlsr  to  those  of  Sokart^^^)  ves  oarrled  out  for  the  oases 
of  iodine  and  earbcm  dioxide  diffusing  throufh  boundary  layer.  Ihe  e<fA- 
atlons  yielded  velocity  and  ooneentration  profiles^  as  well  as  relationships 
between  the  drag  coefficient  and  wall  concentration  or  Inlet  velocity.  Finally^ 
for  each  of  the  wall  conoentrations  speclfLed^  accurate  initial  conditions 
were  obtained  ihlch  would  undoubtedly  serve  as  a  good  Junplng-off  place  for 
a  oaovergence  sdiene  If  the  conplex  eolations  are  solved. 

Once  having  obtained  the  velocity  and  concentration  profiles,  we  must 
not  forget  that  the  original  purpose  of  the  investigation  was  to  exwnlne 
the  stability  character  of  diffusion  boundary  layers.  Here  again,  a  slnpli* 
fled  approach  will  be  used.  It  was  Shown  In  the  TolLnlen-Raylel^  analysis 
that  the  curvature  of  the  velocity  profile  was  fundamental  in  Indicating  the 
stability  of  an  InecDqnressible  flow.  If  the  flow  were  in  the  nature  of  a 
Blaslus  profile,  l.e..  If  It  did  not  exhibit  a  point  of  Inflection,  the  pro* 
file  would  be  stable  regardless  of  Reynolds  mimber.  If  a  profile  showed  a 
definite  point  of  ln£Lecti<m,  as  in  the  ease  of  How  against  an  adverse 
pressure  gradient,  the  boundary  layer  would  be  unstable  for  any  Reynolds 
number.  Although  this  analysis  bolds  only  for  the  case  of  invlscld  flow, 
it  nevertheless  indicates  the  importance  of  the  velocity  profile  In  deter¬ 
mining,  qualitatively  at  least,  the  stability  characteristics  of  the  flow. 

(2) 

Van  Driest  has  extended  the  analysis  to  the  ease  of  eonpresslble  fluids.'  ' 

He  showed  that  if  the  quantity 

(d/dy)  j’p  (du/dy)1  .  0  c^u^  >  1  -  l/^„  (GO) 


-3k. 


a  neutral  or  self-excited  disturbance  appears  in  the  boundary  layer  at  that 

point,  in  th.  dlntrltuUon  01  the  7  Vi  [t/n.]  Is  do- 

eisive  in  determining  the  stabili-^  of  the  boundary  layer.  Ihis  may  be 

seen  by  examining  odn/dy  profiles  for  the  cases  of  heating  and  cooling  the 
(2) 

boundary  layer.'  Ohe  cooling  curves  show  a  tendency  tovard  leveling  out, 
vhereas  the  heating  curves  all  e^tihibit  the  characteristic  hungp,  l.e*,  the 
point  of  inflection.  !Qils,  then,  will  be  one  vay  to  •vmi-tiw  the  stability 
character  of  binary  boundary  layers. 

Another  vay  to  characterize  stability  of  flow  is  to  indicate  the  mini- 
num  critical  Reynolds  number  of  a  system.  !Dils  is  the  smallest  Reynolds 
number  at  which  neutral  disturbances  can  exist.  It  is  obvious  that  we  vould 
like  to  make  the  minimum  critical  Reynolds  number  as  large  as  possible.  One 
way  to  obtain  this  Reynolds  number  is  to  observe  a  plot  of  o  (disturbance 
wavelength)  against  Reynolds  number.  Ihe  curve  Is  usually  paraboloid  In 
shape  and  contains  a  region  in  which  all  self-excited  disturbances  are  aapU- 
fled.  Outside  this  curve  all  disturbances  are  dainped.  Ihe  tip  of  the  para¬ 
boloid  will  be  the  minimum  critical  Reynolds  number.  The  calculation  of 
such  a  curve  Is  very  tedious,  and  Lin^^^^  has  suggested  an  approximation 
fonmila  Wileh  has  been  shown  to  be  within  20  per  cent  of  the  values  obtained 
with  the  more  aoourate  analyses.  The  fOnnula  Is 


where 


(81) 


(82) 
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Lln's  formula  vlU  facilitate  a  qualitative  congparlson  of  the  stability  of 
the  three  systans* 
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vm.  RBSOLIS 

nie  three  systens  investigated  vere  distinguished  by  the  fbreigh  sub- 
stance  used  In  each:  hydrogen^  Iodine  vapor,  or  carbon  dioxide.  Hydrogen 
vas  chosen  because  It  has  been  used  In  blnary-boundaxy-layer.  analyses  previ¬ 
ously,  mostly  on  the  basis  of  Its  ability  ae  a  coolant.  In  this  vay  It  would 
act  as  a  stabilizing  factor.  Howsver,  Its  U^t  weight  would  causa  It  to 
act  as  a  buoying  agent,  tending  to  "lift*  the  boundary  layer  from  Ihe  sur¬ 
face.  Needless  to  say,  this  effect  vould  be  destabilizing  since  It  Increases 
the  boundary-layer  thickness.  Ohe  coolant  effect  of  tiie  hydrogen  vlU  not 
show  up  in  the  analysis,  and  therefore  It  ml^t  be  ezpected  that  Its  sta¬ 
bility  diaraeter  will  not  be  as  good  as  that  of  the  hescvler  materials— 
carbon  dioxide  and  iodine.  Ihese  two  were  diosen  on  -ttie  basis  of  the  Van 
Driest  criteria.  It  was  hoped  that  by  inaki,ng  the  quantity  pda/6y  as  large 
as  possible  near  the  surface,  the  unstable  hump  ole^t  possibly  be  acrolded. 
Also,  the  Introduction  of  heavy  molecules  might  be  thou^t  of  as  "cooling" 
the  boundary  layer  simply  by  Increasing  the  density  at  the  region  near  the 
wall. 

The  equations  were  solved  as  Indicated,  and  the  numerical  results  are 
presented  in  the  Appendices.  An  exmalnatlcn  of  the  velocity  and  concentration 
profiles  bears  out  the  theory  that  hydrogen  acts  as  a  buoying  agent  while 
carbon  dioxide  and  iodine  act  In  the  opposite  direction.  The  hydrogen  dif¬ 
fusing  rapidly  throu^  the  boundary  layer  slows  up  the  a^dlrected  velocity 
of  the  alrstream.  In  Fig.  1,  for  a  wall  concentration  of  hydrogen  of  0.2, 
only  95  per  cent  of  the  free- stream  velocity  has  been  reached  at  11  *  3. 

For  the  same  conditions,  the  free-stream  velocity  for  carbon  dioxide 
and  iodine  injection  is  reached  at  11  =  2.7  and  at  11  =  2.8,  respectively 


Velocity  profiles :  hydrogen 
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(see  Figs*  2  and  3).  Sxs  ame  diaraeterlstlc  of  hydrogen  Is  home  out  hy 
the  concentration  profiles*  Hydrogen  concentrations  extend  out  Into  t2ie 
boundary  layer  well  past  T)  ■  3  (Fig*  4),  vhereas  the  concentration  of  carbon 
dioxide  and  Iodine  for  ■  3  are  fractlcally  zero  (Figs*  ^  and  6)* 

Eolation  (53)  gives  the  relation  between  the  Injection  rate  and  the 
concentration  of  the  foreign  material*  Because  of  Its  hl^  diffusion  coe> 
fflclent  snd  low  density^  hydrogen  should  have  a  hlc^  velocity  and  a  lov 
Injection  rate*  Examination  of  Fig*  7  Indicates  that  this  Is  generally 
true*  OBie  Injection  rates  of  the  three  materials  are  alioost  the  same  at  a 
vail  concentration  of  0*2;  but  It  must  be  remembered  that  even  at  this  point 
the  presence  of  the  forel^i  constituents  Is  beginning  to  taper  off,  and  for 
smaller  eoneentratlona  ve  may  expect  the  Injection  rate  to  drop  to  zero 
quickly*  At  the  hlgh«r  concentrations,  however,  it  Is  seen  that  large  diffai> 
ences  occur  in  the  Injection  rates  of  the  three  materials*  !Ihe  essential 
nature  of  the  curves  of  Fig*  7  !•  similar  to  a  plot  of  the  density  ratios 
of  the  three  materials  (Fig*  8)*  Oharefore,  considering  that  the  concent 
tratlon  gradients  at  the  surfhce  are  almost  the  saae  for  all  materials.  It 
Is  seen  that  the  Injection  rate  Is  dependant  mainly  upon  the  density  ratio* 

At  lov  oonoentratlons  all  three  density  ratios  apporoadi  unity,  vhereupon  the 
Schmidt  xunber  or  diffusion  coefficient  becomes  the  Important  factor*  At 
extremely  lov  vail  cencen'brations  the  Injection  rates  at  the  vail  become 
almost  proportional  to  the  diffusion  coefficients. 

Ohe  curves  in  Fig*  9  shoving  the  relationship  between  the  drag  coef¬ 
ficient  and  the  Injection  rata  (vhare  f^  Is  a  maasure  of  Injection  rate,  as 
ahovn  by  Eq.  (52))  indicate  that  hydrogen  is  superior  to  carbon  dioxide 
or  Iodine  as  a  fidetlon-reduoing  agent*  If  separation  ocoun  vhen  the 
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ocity  profiles;  carbon  dioxide 


Fig. 3  —  Velocity  profiles:  iodine 
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P* 


Fig. 4  — Concentration  profiles-,  hydrogen 
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Fig.  5  —  Concentrotion  profiles:  carbon  dioxide 


,43 


Fig. 6  —  Concentrotion  profiles:  iodine 


Fig.  7 —  Injection  rot 


Fig.  8  —  P  vs  C| 
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Fig.  9 —  Friction  coefficient  vs  injection  rote 


-47- 


firletlon  coefficient  is  zero^  extending  -ttie  hydrogen  curve  shovs  that  at 
f^  a  0.3$  the  hydrogen  would  have  blown  the  boundary  layer  off  the  surface; 

In  the  case  of  carbon  dioxide  or  iodine  it  la  still  firmly  entrenched.  This 
indicates  the  destabilizing  character  of  hydrogen  injection.  Both  carbon 
dioxide  and  iodine  increase  the  friction  coefficient  over  that  value  ihlch 
it  would  have  with  air  injection.  At  lov  concentrations  this  increase  is 
seen  to  be  small.  Carbon  dioxide,  particularly,  gives  values  which  differ 
by  a  small,  aliaost  constant  factor  froa  the  normal  air  values.  In  the  case 
of  Iodine,  the  breach  becomes  geometrically  greater  as  the  Injection  rate 
increases,  ^e  drag  coefficient  is  plotted  against  wall  concentration  in 
Fig.  10. 

Ihe  next  series  of  graphs  (Figs.  11  -  13)  shova  the  "stability"  functloin, 
p  d/d|^u/u^~|  plotted  against  T|  fbr  varying  wall  concentrations.  As  was 
stated  previously,  the  characteristic  huiiqp  indicative  of  a  point  of  in¬ 
flection  shows  that  the  systam  Is  unstable.  On  the  other  hand.  It  is  to 
be  noted  that  this  criterion  applies  to  an  incoapiresslble  fluid  and  that 
ours  Is  not  truly  incoBqpresslble,  since  the  density  and  viscosity  vary 
throughout  the  boundary  layer.  The  variation  of  these  two  parameters  Is 
dependent  upon  concentration  and  xx>t  upon  tmaparature.  Ihexefore,  any 
conclusions  ^ich  can  be  drawn  are  at  beat  <iualitative  and  are  Intended  to 
show  direction  rather  than  order  of  magnitude.  Since  for  large  enough  i]  all 
the  stability  curves  must  approach  zero,  a  hump  will  automatically  appear 
if  the  slbpe  of  the  function  at  the  wall  is  greater  than  zero.  In  the 
case  of  hydrogen  it  is  lasedlately  spparent  that  instability  exiata  for 
all  oonditions.  Itr  oonoentrations  of  0.2  and  0.4  the  bmps  are  still 
(^Varent,  but  far  hl^er  oonoentrations  they  appear  In  the  boundary  layer 
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Fig.  10 — Friction  coefficient  vs  injection  composition 
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Fig.  II —  P  "uT  vs  “J?:  hydrogen 


Fig.  12 —  P  "uT  vs  ^7;  carbon  dioxide 


iodine 
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fOir  T)  >  3*  In  case  of  carbon  dloadda.  It  can  be  seen  that  the  infiee- 
tiOD  points  have  moved  tovard  the  left;  for  a  ocnoentratlon  of  O.03  the  Idp 
flection  husqp  has  almost  disappeared,  m  the  ease  of  Iodine,  the  Inflection 
points  are  close  to  the  left  axis  and  are  no  Ioniser  discernible  for  the 
three  lover  conoentratlans.  Even  for  a  eenoentration  of  0.2,  calculations 
shov  that  stability  has  bean  achieved.  lOr  lover  concentration  It  is  evi¬ 
dent  that  these  curves  hacve  flattened  out  and  that  they  shov  a  stable  charac¬ 
ter.  A  careful  numerical  analysis  should  indicate  the  most  advantageous 
concentration.  All  three  materials  tend  toward  the  Blaslus  preflle  In  the 
limit  (as  concentration  approadies  zero).  At  large  concentrations  aU  are 
obviously  unstable.  It  vlU  be  In  the  range  of  small  Injection  rates  that 
the  difference  In  the  three  materials  vUl  be  most  noticeable.  Oils  vould 
actually  be  the  ease  If  an  eaqperlmental  model  vere  to  be  made.  It  goes 
vlthout  saying  that  a  careful  numerical  analysis  of  the  lover  concentrations 
Is  In  order. 

Oie  critical  Reynolds  numbers  are  plotted  in  Fig.  lU.  Here  it  Is  again 
evident  that  hydrogen  injection  reduces  the  stability  character  of  a  bounda¬ 
ry  layer  significantly,  even  at  low  concentrations.  Iodine,  of  the  three 
materials,  gives  the  leut  reduction  azid  Is  the  best  ixijection  medium  sic- 
cording  to  stability  considerations.  A  general  direction  of  comparative 
stability  characteristics  has  been  shovn. 

Injection  media,  such  as  iodine,  that  are  less  destabilizing  prevent 
hmediate  transition  and  Increase  the  area  affected  by  the  lovh.drag  laminar 
flow.  Hovever,  the  medium  Itself  produces  an  increase  of  dreg  (above  that 
obtained  xuader  similar  conditions  vlth  hydrogen  Injection)  so  that  some  of 
Its  advantage  is  negated.  Ihe  vail  concentration  at  vhlch  iodine  most 


Fig.  14  —  The  influence  of  moss  transfer  on  the  stability 
of  the  incompressible  laminar  boundary  layer 


effiectlvely  prolonga  laminar  flow  must  first  be  obtained.  Qie  critical 
Reynolds  number  Hor  this  condition  may  then  be  found,  and  a  conqparatlve 
estimate  may  be  made  between  turbulent,  IsQWimeable-llat-plate,  and  Iodine- 
air  boundary  layers* 
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Appendlx  A 

SOLUTION  OF  OHE  LAMINAR  BOUNDARY  LAYER 
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Tadile  2 

SOLUIlOH  OF  THE  LAMINAR  BOUNDARY  LAYER  WITH 
CARB(»I  DlOKinS  INJECTION 
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(Ci„  -  0.001) 

0 

0.2 

0.4 

0.6 

0.8 

— 

0 

0.1298 

0.2589 

0.3B54 

0.5062 

0 

0.1327 

0.2646 

0.3936 

0.5166 

-0.0007 

0.0259 

0.1054 

0.2372 

0.4195 

1.3273 

1.3251 

1.3090 

1.2660 

1.1865 

10,0 

8,7018 

7.4110 

6.1464 

4.9378 

-6,4950 

-6.4851 

-6.4109 

-6.2128 

-5.8454 

1.0 

1.2 

1.4 

1.6 

1.8 

0.6179 

0.7167 

0.7998 

0.8^ 

0.9150 

0.6296 

0.7288 

0.8114 

0.8760 

0.9232 

0.6492 

0.9214 

1.2301 

1.5682 

1.9286 

1.0670 

0.9125 

0.7363 

0,5568 

0.3926 

3.8212 

2.8332 

2.0020 

1.3419 

0.8498 

-5.2905 

-4.5673 

-3.7330 

-2,8708 

-2.0672 

2.0 

2.2 

2.4 

2.6 

2.8 

2.9 

0.9493 

0.9716 

0.9851 

0.9927 

0.9^ 

0.9978 

0.9555 

0.9758 

0.9878 

0.9942 

0.9975 

0.9^ 

2.30W 

2.6914 

3.0^ 

3.1i809 

3.8793 

4.0789 

0.2572 

0.1560 

0.0876 

0.0454 

0.0218 

0.0146 

0.5066 

0.2836 

0.1488 

0.0730 

0.0334 

0.0221 

-I.3B88 

-0.8684 

-0.5046 

-0.2723 

-0.1363 

-0.0938 

-56. 


Table  2  (Contd.) 


T1 

°iw-°i 

u 

0; 

yo 

^2 

X  10^ 

X  io3 

(Ci^  -  0.005) 

0 

-0.0035 

1.3236 

5.0 

-3.2411 

0.2 

0.1324 

0.0231 

1.3222 

4.3529 

-3.2378 

o.u 

0.2641 

0.1025 

1.3068 

3.7089 

-3.2025 

0.6 

0.3930 

0.2342 

1.2646 

3.0776 

-3.1054 

0.8 

0.5052 

0.5159 

0.4164 

1.1859 

2.4738 

-2.9235 

1.0 

0.6169 

0.6289 

0.6460 

1.0672 

1.9155 

-2.6477 

1.2 

0.7158 

0.7282 

0.9180 

0.9132 

1.4210 

-2.2874 

l.U 

0.7990 

0.8108 

1.2265 

0.7373 

1.0043 

-1.8709 

1.6 

0.8652 

0.8755 

1.5645 

0.5580 

0.6739 

-1.2327 

1.8 

0.9146 

0.9229 

1.9248 

0.3937 

0.4270 

-1.0376 

2.0 

0.9490 

0.9552 

2.3010 

0.2581 

0.2548 

-0.6976 

2.2 

0.9714 

0.9757 

2.6875 

0.1567 

0.1427 

-0.4366 

2.4 

0.9850 

0.9877 

3.0804 

0.0880 

0.0749 

-0.2539 

2.6 

0.9926 

0.9942 

3.4770 

0.0457 

0.0368 

-0.1371 

2.8 

0.9^ 

0.9974 

3.8754 

0.0219 

0.0169 

-0.0687 

2.9 

0.9978 

0.9984 

4.0750 

0.0147 

0.0111 

-0.0473 

■  0.01) 


0 

[ - 

0 

0 

-0.0069 

1.3190 

10.0 

-6.4661 

0.2 

0.1289 

0.1321 

0.0195 

1.3184 

8.7106 

-6.4638 

0.4 

0.2574 

0.2634 

0.0989 

1.3041 

7.4264 

-6.3976 

0.6 

0.3834 

0.3922 

0.2305 

1.2629 

6.1663 

-6.2060 

0.8 

0.5040 

0.5150 

0.4125 

1.1851 

4.9598 

-5.8491 

1.0 

0.6157 

0.6280 

0.6419 

1.0673 

3.8433 

-5.3017 

1.2 

0.7147 

0.7273 

0.9137 

0.9141 

2.8534 

-4.5840 

1.4 

0.7981 

0.8101 

1.2221 

0.7387 

2.0191 

-3.7528 

1.6 

0.8645 

0.8750 

1.5599 

0.5595 

1.3554 

-2.8908 

1.8 

0.9140 

0.9225 

1.9201 

0.3952 

0.8596 

-2.0852 

2.0 

0.9487 

0.9550 

2.2961 

0.2593 

0.5132 

-1.4033 

2.2 

0.9712 

0.9755 

2.6826 

0.1576 

0.2878 

-0.8790 

2.4 

0.9^9 

0.9676 

3.0755 

0.0886 

0.1512 

-0.5116 

2.6 

0.9926 

0.9941 

3.4720 

0.0460 

0.0743 

-0.2765 

2.8 

0.9^ 

0.9974 

3.8704 

0.0221 

0.0341 

-0.1387 

2.9 

0.9977 

0.9983 

4.0700 

0.0148 

0.0225 

-0.0955 
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Table  2  (Contd.) 


T\ 

Ciw  ■  Cl 
ClW 

u 

Ce 

yo 

CM 

^3 

y4 

II 

> 

#-i 

0 

0.05) 

0 

0 

0 

-0.0355 

1.2810 

0.050 

-0.031670 

0.2 

0.1250 

0.1292 

-0.0092 

1.2879 

0.043749 

0.4 

0.2504 

0.25^ 

0.0693 

1.2813 

0.037479 

-0.031676 

0.6 

0.3744 

0.3856 

0.1998 

1.2482 

0.031281 

-0.030918 

0.8 

0.4940 

0.5076 

0.3^4 

1.1786 

0.025301 

-0.029314 

1.0 

0.6056 

0.6206 

0.6081 

1.0682 

0.019721 

-0.026751 

1.2 

0.7053 

0.7205 

0.8784 

0.9210 

0.014733 

-0.023296 

1.4 

0.7901 

0.8042 

1.1852 

0.7495 

0.010493 

-0.019214 

1.6 

0.8582 

0.8704 

1.5218 

0.5719 

0.007090 

-0.014914 

1.8 

0.9094 

0.9192 

1.8809 

0.4070 

0.004527 

-0.010842 

2.0 

0.9456 

0.9527 

2.2562 

0.2691 

0.002722 

-0.007554 

2.2 

0.9693 

0.9741 

2.6422 

0.1649 

0.001537 

-0.004^2 

2.4 

0.9837 

0.9868 

3.0348 

0.0935 

0.000813 

-0.002723 

2.6 

0.9920 

0.9937 

3.4311 

0.0490 

0.000402 

-0.001483 

2.8 

0.9963 

0.9972 

3.8294 

0.0237 

0.000186 

-0.000750 

2.9 

0.9975 

0.9982 

4.0290 

0.0160 

0.000123 

-0.000518 

^Ciw  “ 

0.2) 

0 

0 

0 

-0.1539 

1.1264 

0.2000 

-0.1157 

0.2 

0.1097 

0.1168 

-0.1291 

1.1604 

0.1780 

-0.1187 

0.4 

0.2226 

0.2363 

-0.0543 

1.1633 

0.1554 

-0.1208 

0.6 

0.3378 

0.3567 

0.0704 

1.1824 

0.1324 

-0.1207 

0.8 

0.4522 

0.4750 

0.2440 

1.1463 

0.1096 

-0.1174 

1.0 

0.5625 

0.5873 

0.4641 

1.0679 

0.0875 

-0.1102 

1.2 

0.6646 

0.6893 

0.7267 

0.9477 

0.0671 

-0.0968 

1.4 

0.75^^5 

0.7772 

1.0264 

0.7949 

0.0491 

-0.0841 

1.6 

0.8293 

0.8486 

1.3568 

0.6260 

0.0342 

-0.0675 

1.8 

0.8878 

0.9030 

1.7110 

0.4603 

0.0224 

-0.0507 

2.0 

0.9305 

0.9416 

2.0827 

0.3148 

0.0139 

-0.0356 

2.2 

0.9596 

0.9670 

2.4662 

0.1997 

0.0081 

-0.0232 

2.4 

0.9780 

0.9826 

2.8573 

0.1172 

0.0044 

-0.0141 

2.6 

0.9^ 

0.9915 

3.2527 

0.0636 

0.0022 

-0.0080 

2.8 

0.9946 

0.9961 

3.6506 

0.0319 

0.0011 

-0.0042 

2.9 

0.9^ 

0.9974 

3.8500 

0.0219 

0.0007 

-0.0029 
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Table  2  (Contd.) 


T1 

■^1 

U 

^0 

^2 

^3 

^4 

0.4) 

0 

0 

0 

-0.3474 

0.8850 

0.4000 

0.2 

0.0876 

0.0959 

-0.3257 

0.9502 

0.3649 

0.4 

0.1813 

0.1979 

-0.2593 

1.0106 

0.3275 

-0.2176 

0.6 

0.2806 

0.3050 

-0.1472 

1.0549 

0.2878 

-0.2253 

0.8 

0.3838 

0.4148 

0.0112 

1.0706 

0.2465 

-0.223B 

1.0 

0.4884 

O.523B 

0.2146 

1.0470 

0.2046 

-0.2238 

1.2 

0.5906 

0.6275 

0.4605 

0.9779 

0.1637 

-0.2110 

1.4 

0.6861 

0.7215 

0.7446 

0.8658 

0.1255 

-0.1895 

1.6 

0.7708 

0.8020 

1.0615 

0.7216 

0.0917 

-0.1610 

1.8 

0.8413 

0.8668 

1.4048 

0.5629 

0.0635 

-0.1285 

2.0 

0.8963 

0.9156 

1.7682 

0.4092 

0.0415 

-0.0959 

2.2 

0.9363 

0.9496 

2.1460 

0.2764 

0.0255 

-0.0666 

2.4 

0.9633 

0.9720 

2.5333 

0.1729 

0.0147 

-0.0431 

2.6 

0.9802 

0.965^^ 

2.9265 

0.1001 

0.0079 

-0.02^ 

2.8 

0.9900 

0.9930 

3.3231 

0.0536 

0.0040 

-0,0144 

2.9 

0.9931 

3.5222 

0.03^ 

0.0028 

-0.0104 

^‘^Iw  •  ' 

0.6) 

0 

0 

0 

-0.5990 

0.5955 

0.6000 

-0.2252 

0.2 

0.0630 

0.0681 

-0.5824 

0.6762 

0.5622 

-0.2563 

0.4 

0.1329 

0.1447 

-0.5306 

0.7619 

0.^202 

-0.2683 

0.6 

0.2102 

0.2299 

-0.4408 

0,8452 

0.4739 

-0.2896 

0.8 

0.2946 

0,3228 

-0.3106 

0.9154 

0.4232 

-0.3076 

1.0 

0.3852 

0.4212 

-0,1388 

0.9602 

0,3689 

-0.3192 

1.2 

0.4790 

0.5218 

0.0745 

0.9674 

0.3669 

-0.3192 

1.4 

0.5752 

0.3275 

0.9287 

0.2549 

-0.3105 

1.6 

0.6673 

0.6164 

0.8435 

0.1996 

-0.2862 

1.8 

0.7516 

0.9362 

0.7204 

0.1491 

-0.2496 

2.0 

0.8241 

0.8571 

1.2813 

0.5756 

0.1055 

-0.2045 

2.2 

0.8825 

0.9078 

1.6456 

0.4265 

0.0705 

-0.1566 

2.4 

0.9262 

0.9440 

2.0237 

0.2962 

0.0443 

-o.m7 

2.6 

0.9565 

0.9682 

2.4m 

0.1698 

0.0261 

-0.0740 

2.8 

0.9760 

0.9631 

2.8043 

0.1126 

0.0144 

-0.0455 

2.9 

0.9826 

0.9680 

3.0022 

0.0842 

0.0104 

-0.03>^ 

0 

0 

-0.9393 

0.8000 

0.031^5 

0.0319 

-0.9309 

0.7724 

0.0742 

0.0709 

-0.9036 

0.7406 

0.1201 

0.1182 

-0.8538 

0.4753 

0.7039 

0.1728 

0.1748 

-0.7779 

0.5678 

0.6618 

0.2329 

0.2412 

-0.6720 

0.6644 

0.6136 

0.3009 

0.3175 

-0.5327 

0.7565 

0.5593 

0.3763 

0.4022 

-0.3574 

0.8322 

0.4990 

0.4580 

0.4928 

-0.1448 

0.8781 

0.4336 

0.5440 

0.5854 

0.1048 

0.8824 

0.361(8 

0.6306 

0.6752 

0.3890 

0.8388 

0.2955 

0.7139 

0.7573 

0.7039 

0.7499 

0.2289 

0.7893 

0.8280 

1.0444 

0.6275 

0.1685 

0.8533 

0.8847 

1.4048 

0.48^ 

0.1173 

0.9039 

0.9273 

1.7800 

0.3548 

0.0769 

0.9240 

0.9436 

1.9716 

0.2937 

0.0608 

-0.1766 

-0.2009 

-0.2283 

-0.2587 

-0.2911 

-0.3239 

-0.35^ 

-0.3^1 

-0.3957 

-0.3971 

-0.^12 

-0.3473 

-0.2980 

-0.239*^ 

-0.1792 

-0.1509 
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Table  3 

SOLUnOI  OF  TBl  LAMUAR  BOUHOARX  LAXER  WITH 
HYDROOEH  IMJECTIOH 


11 

■M 

U 

yo 

_ 

^2 

73X10*^ 

y^xlO*^ 

-  0.0001) 

0 

0 

0 

-0.0002 

1.3278 

1.0 

-0.39761 

0.2 

0.0796 

0.1327 

0.0263 

1.3250 

0.92088 

-0.39748 

0.4 

0.1591 

0.2645 

0.1058 

1.3088 

0.84087 

-0.39639 

0.6 

0.2382 

0.3933 

0.2314 

1.2657 

0.76176 

-0.39347 

0.8 

0.3165 

0.5165 

0.41^ 

I.IB62 

0.68350 

-0.36786 

1.0 

0.3933 

0.6295 

0.6488 

1.0669 

0.60670 

-0.37887 

1.2 

0.4679 

0.7287 

0.9208 

0.9125 

0.53209 

0.46047 

-0.36600 

1.4 

0.5395 

0.6074 

0.8112 

1.2291 

0.7364 

-0.34904 

1.6 

0.8758 

1.5670 

0.5570 

0.39265 

-0.32815 

1.8 

0.6706 

0.9231 

1.9271 

0.3929 

0.32937 

-0.30377 

2.0 

0.7287 

0.9554 

2.3031 

0.2574 

0.27127 

-0.27666 

2.2 

0.7812 

0.9758 

2.6896 

0.1562 

0.21879 

0.24776 

-0.2Ij^ 

2.4 

0.8278 

0.9877 

3.0824 

0.0877 

0.17219 

2.6 

0.8685 

0.99^ 

3.4788 

0.0455 

0.13152 

-0.18864 

2.8 

0.9034 

0.9974 

3.8772 

0.0218 

0.096^ 

-0.16034 

2.9 

O.91B7 

0.9984 

4.0768 

0.0146 

0.08128 

-0.145^ 

-  0.001) 


0 

0 

0 

-0.0017 

1.3196 

10.0 

-3.8119 

0.2 

0.0772 

0.1320 

i  0.0243 

1.3178 

9.2278 

-3.8107 

0.4 

0.1543 

0.2632 

0.1025 

1.3023 

8.4573 

-3.8007 

0.6 

0.2309 

0.3916 

0.2322 

1.2605 

7.6913 

-3.7731 

0.8 

0.3066 

0.5141 

0.4118 

1.1826 

6.9341 

-3.7199 

1.0 

0.3808 

0.6269 

0.6384 

1.0657 

6.1916 

-3.6345 

1.2 

0.4529 

0.7261 

0.9074 

0.9136 

5.4706 

-3.5120 

1.4 

0.5221 

0.8089 

1.2128 

0.7395 

4.7790 

-3.3505 

1.6 

0.5876 

0.8739 

1.5480 

0.5613 

4.1243 

-3.1514 

1.8 

0.6486 

0.9216 

1.9058 

0.3975 

3.5136 

-2.9188 

2.0 

0.7047 

0.9543 

2.2796 

0.2616 

2.9529 

-2.6600 

2.2 

0.7553 

0.9751 

2.6646 

0.1596 

2.4465 

-2.2422 

2.4 

0.8003 

0.9873 

3.0562 

0.0900 

1.9967 

-2.0999 

2.6 

0.8396 

0.9940 

3.4516 

0.0470 

1.6041 

2.8 

0.8733 

0.9978 

3.8492 

0.0226 

1.2672 

-1.5467 

2.9 

0.8881 

0.9983 

4.0484 

0.0152 

1.1188 

-1.4172 
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Table  3  (Contd.) 


T1 

B!!B9 

yo 

-  -  1 

^2 

^3 

X  10^ 

0.05) 

0 

0 

0 

-0.0631 

1.0903 

0.05 

-1.3245 

0.2 

0.0874 

0.1134 

-0.0492 

1.1036 

0.045632 

-1.3301 

O.U 

0.1719 

0.2273 

-0.0062 

1.1105 

0.041404 

-1.3330 

0.6 

0.2535 

0.3410 

0.0681 

1.1041 

0.037324 

-1.3304 

0.8 

0.3319 

0.4526 

0.1757 

1.0776 

0.033406 

-1.3202 

1.0 

0.4066 

0.5599 

0.3180 

1.0253 

0.029668 

-1.3001 

1.2 

0.4774 

0.6602 

0.4958 

0.9446 

0.026129 

-1.2686 

1.4 

0.5438 

0.7506 

0.7081 

0.8367 

0.022810 

-1.2247 

1.6 

0.6054 

0.8268 

0.9540 

0.7077 

0.019727 

-1.1684 

1.8 

0.6621 

0.8932 

1.2307 

0.5680 

0.016896 

-1.1003 

2.0 

0.7135 

0.9433 

1.5346 

0.4302 

0.014324 

-1.0220 

2.2 

0.7596 

0.98^ 

1.8617 

0.3058 

0.012018 

-0.9357 

2.4 

0.8005 

1.0055 

2.2079 

0,2033 

0.009973 

-0.8440 

2.6 

0.8363 

1.0218 

2.5695 

0.1259 

0.008184 

-0.7497 

2.8 

0.8672 

1.0316 

2.9431 

0.0724 

0.006638 

-0,6558 

2.9 

0.8810 

1.0347 

3.1336 

0.0534 

0,005952 

-0.60^ 

0.2) 

0 

0 

0 

-0.1462 

0.6429 

0.2000 

-2.5853 

0.2 

0.0934 

0.0792 

-0.1418 

0.6658 

0,1813 

-2.6451 

0.4 

O.1B23 

0.1594 

-0.1274 

0,6876 

0.1635 

-2.6963 

0.6 

0.2665 

0.2ll04 

-0.1014 

0.7063 

0.1467 

-2.7413 

0.8 

0.3457 

0.3219 

-0.0622 

0.7198 

0.1309 

-2.7704 

1.0 

0.4197 

0.4030 

-0.0077 

0.7257 

0.1160 

-2.7824 

1.2 

0.4885 

0.4829 

0.0641 

0.7214 

0.1023 

-2.7744 

1.4 

0.5518 

0.5606 

0.1552 

0.7047 

0.0896 

-2.7444 

1.6 

0.6098 

0.6346 

0,2676 

0.6738 

O.OTto 

-2.6909 

1.8 

0.6624 

0.7037 

0.4030 

0.6278 

0.0675 

-2.6137 

2.0 

0.7096 

0.7665 

0.5625 

0.5676 

0.0580 

-2.5133 

2.2 

0.7521 

0.8218 

0.7469 

0.4957 

0.0496 

-2.3915 

2.4 

0.7897 

0.8689 

0.9562 

0.4160 

0.0421 

-2.2508 

2.6 

0.8228 

0.9073 

1.1899 

0.3341 

0.0354 

-2.0^ 

2.8 

0.8517 

0.9373 

1.4467 

0.2554 

0.0296 

-1.9264 

2.9 

0.86i«8 

0.94^ 

1.5832 

0.2190 

0,0270 

-1.8392 

-64- 


Teible  3  (Cantd.) 


n 

mm 

imi 

U 

^0 

^2 

^3 

74  X  10^ 

>  0.4) 

0 

0 

0 

-0.2044 

0.4U6 

0.4000 

-2.7106 

0.2 

0.0860 

0.0618 

-0.2024 

0.4368 

0.3656 

-2.8654 

0.4 

0.1698 

0.1254 

-0.1958 

0.4622 

0.3321 

-3.0134 

0.6 

0.2509 

0.1906 

-0.1838 

0.4870 

0.2996 

-3.1498 

0.8 

0.3285 

0.2571 

-0.1652 

0.5103 

0.2686 

-3.2700 

1.0 

0.4020 

0,3244 

-0.1388 

0.5308 

0.2392 

-3.3693 

1.2 

0.4709 

0.3922 

-0.1033 

0.5473 

0.2116 

-3.4436 

1.4 

0.5347 

0.4597 

-0.0570 

0.5583 

0.1861 

-3.4896 

1.6 

0,593l^ 

0.5263 

0.0020 

0,5621 

0.1626 

-3.5046 

1.8 

0.6467 

0.5913 

0.0754 

0.5573 

0.1413 

-3.  *♦873 

2.0 

0.6948 

0.6536 

0.1653 

0,5425 

0.1221 

-3.4370 

2.2 

0.7378 

0.7124 

0,2735 

0.5169 

0.1049 

-3. 35*^6 

2.4 

0,7760 

0.7668 

0.4018 

0.4805 

0.0896 

-3.2418 

2.6 

0.8096 

0.8159 

0.5516 

0,4340 

0.0762 

-3.1012 

2.8 

0.8390 

0.8590 

0.7239 

0.3795 

0.0644 

-2.9364 

2.9 

0.8522 

0.8781 

0.8186 

0.3500 

0.0591 

-2.8462 

0.6) 

0 

0 

0 

-0.2468 

0.2654 

0.6000 

-2.1815 

0.2 

0.0674 

0.0467 

-0.2457 

0.2864 

0.5596 

-2.4018 

0,4 

0.1367 

0,0960 

-0.2423 

0.3125 

0.5180 

-2.6262 

0.6 

0.2070 

0.1478 

-0.2359 

0.3373 

0.4758 

-2.8493 

0.8 

0.2774 

0.2018 

-0.2261 

0.3623 

0.4336 

-3.0652 

1.0 

0,3467 

0.2578 

-0.2121 

0.3869 

0.3920 

-3.2670 

1.2 

0.4141 

0.3154 

-0.1930 

0.4102 

0.3516 

-3,4484 

1.4 

0.4785 

0.3740 

-0.1678 

0,4315 

0.2769 

-3.7261 

1.6 

0.5393 

0.4334 

-0.1353 

0,4495 

0.2769 

-3.7261 

1.8 

0.5959 

0.4928 

-0.0940 

0.4633 

0.2424 

-3.8130 

2.0 

0.6479 

0.5516 

-0.0425 

0.4714 

0.2113 

-3.8609 

2.2 

0.6950 

0.6092 

0.0211 

0.4728 

0.1830 

-3.8686 

2.4 

0,7374 

0.6648 

0,0986 

0.4663 

0.1576 

-3.8361 

2.6 

0.7751 

0.7177 

O.I9IB 

0.4511 

0.1350 

-3.7647 

2.8 

0.8083 

0.7671 

0.3026 

0.4^ 

0.1150 

-3.6570 

2.9 

0.8234 

0.7902 

0.3652 

0.4115 

0.1060 

-3.5905 
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Tibla  3  (Coatd.) 


(C^y  -  0.8) 


0 

0 

0 

-0.2871 

0.1480 

0.8000 

-1.2689 

0.2 

0.0396 

0.0296 

-0.2866 

0.1650 

0.7683 

-1.4698 

O.U 

0.0836 

0.0621 

-0.2859 

0.1836 

0.7332 

-1.69^ 

0.6 

0.1317 

0.0974 

-0.2819 

0.2036 

0.6946 

-1.9334 

0.8 

O.IS37 

0.1357 

-0.2773 

0.6530 

-2.1906 

1.0 

0.2390 

0.1769 

-0.2708 

0.2476 

0.6088 

-2.4585 

1.2 

0.2967 

0.2208 

-O.261B 

0.2710 

0.5626 

-2.7306 

1.4 

0.3560 

0.2672 

-0.2500 

0.2^ 

-2.9990 

1.6 

0.4157 

0.3159 

-0.2347 

0.^84 

-3.2550 

1.8 

0.4746 

0.3664 

-0.2152 

0.3412 

0.4203 

-3.4899 

2.0 

0.531fl 

0.4183 

-0.1905 

0.3^ 

0.3746 

-3.69^ 

2.2 

0.5862 

0.4711 

-0.1595 

O.3B08 

0.3310 

-3.80^ 

2.4 

0.6372 

0.5242 

•0.12U 

0.3958 

0.2902 

-3.9920 

2.6 

0.6843 

0.5770 

-0.0738 

0.4062 

0.2526 

-4.0740 

2.8 

0.7270 

0.6290 

-0.0159 

0.4110 

0.2184 

-4.1089 

2.9 

0.7467 

0.6544 

0.0176 

0.4109 

0.2026 

-4.1088 
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Table 


80LOTI0I  OP  IHl  LAKHAR  BIIARI  BOOTDARI 
.  um  TAnm  TB.'nEnvm 


0 

0.2 

0.1^ 

0.6 

0.8 


1.6 

1.8 


0 

0.2 

0.4 

0.6 

0.8 


0 

0.1573 

0.3130 

0.4633 

0.6022 

0.7232 

0.8212 

O.893B 

0.9426 

0.9720 

0.9677 

0.9952 

0.9983 

0.9995 

0.9996 

0.9999 


1.8 

2.0 

2.2 

2.4 

2.6 

2.8 

2.9 


0 

0.1572 

0.3128 

0.4630 

0.6020 

0.7230 

0.8210 

0.8937 

0.9425 

0.9719 

0.9877 

0.9952 

0.9963 

0.9995 

0.9998 

0.9999 


u 

^9 

y  X  lo'* 

% 

*0 

d 

74 


X  10^ 


(®1 


0 

0.1328 

0.2647 

0.3936 

0.5167 

0.6296 

0.72^ 

0.8115 

0.8761 

0.9233 


0 

0.0265 

0.1061 

0.2379 

0.4203 

0.6500 

0.9223 

1.2310 

1.5691 

1.9295 

2.3057 

2.6923 

3.0853 

3.4818 

3.8803 

4.0799 


1.3282 

1.3258 

1.3095 

1.2664 

I.1B67 

1.0670 

0.9124 

0.7360 

0.5565 

0.3924 

0.2569 

0.1559 

O.C875 

0.0454 

0.0217 

0.0146 


1.0 

0.84269 

0.68695 

0.53670 

0.39761 

0.27675 

0.17679 

0.10616 

0.05740 

0.02804 

0.01230 

0.00482 

0.00168 

0.00053 

0.00015 

0.00008 


0 

0.1327 

0.2646 

0.3937 

0.5167 

0.6297 

0.7289 

0.8114 

0.8760 

0.9233 

0.9555 

0.9759 

0.9878 

0.9942 

0.9^ 

0.9969 


-0.0005 
0.0261 
0.1057 
0.2375 

0.4200 

0.6497 

0.9220 

1.2307 

1.5688 

1.9292 

2.3054 

2.6920 

3.0850 

3.4816 

3.6806 

4.0796 


1.3278 

1.3256 

1.3094 

1.2663 

I.IB67 

1.0671 

0.9125 

0.7362 

0.5567 

0.3925 

0.2570 

0.1560 

0.0875 

0.0454 

0.0317 

0.0146 


10.0 
8.4281 
6.8716 

5.3^ 

3.9605 

2.7696 

1.7694 

1.0626 

0.5746 

0.2808 

0.1232 

0.0483 

0.0169 

0.0053 

0.0028 

0.0008 


-0.78718 

-0.78492 

-0.76925 

-0.72833 

-0.65511 

-0.55092 

-0.1^683 

-0.33076 

-O.190& 

-0.10786 

-0.05410 

-0.02396 

-0.00934 

-0.00320 

-0.00096 

-0.00050 


.7.8708 

-7.8494 

-7.6937 

-7.2856 

-6.5542 

-5.5126 

-4.2716 

-3.0103 

-1.9088 

-1.0796 

-0.5417 

-0.2399 

-0.0936 

-0.0321 

-0.0097 

-0.0050 
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Tal>le  U  (Contd. ) 


n 

■Wif 

u  ! 

X  10^  1 

X  10^ 

■a 

■ai 

^0 

^2 

0.005) 

0 

0 

0 

-0.0024 

1.3260 

5.0 

-3.9332 

0.2 

0.1566 

0.1325 

0.0242 

1.3243 

4.2168 

-3.9250 

0.4 

0.3119 

0.2643 

0.1039 

1.3086 

3.4403 

-3.6496 

0.6 

0.4620 

0.3933 

0.2359 

1.2660 

2.6901 

-3.6479 

0.8 

0.6009 

0.5163 

0.4185 

1.1868 

1.9956 

-3.2839 

1.0 

0.7221 

0.6293 

0.6483 

1.6675 

1.3894 

-2.7638 

1.2 

0.8203 

0.72^ 

0.9207 

0.9131 

0.8983 

-2.1430 

1.4 

0.8932 

0.8112 

1.2294 

0.7369 

0.5337 

-1.5111 

1.6 

0.9422 

0.8758 

1.5675 

0.5573 

0.2888 

-0.9587 

1.8 

0.9718 

0.9232 

I.92&O 

0.3920 

0.1412 

-0.5426 

2.0 

0.9676 

0.9554 

2.3042 

0.2575 

0.0620 

-0.2723 

2.2 

0.9951 

0.9758 

2.6908 

0.1563 

0.0243 

-0.1207 

2.4 

0.9983 

0.9876 

3.0837 

0.0877 

0.0085 

-0,0471 

2.6 

0.9995 

0.9942 

3.4802 

0.0455 

0,0026 

-0.0162 

2.8 

0.9998 

0.9975 

3.8786 

0.0218 

0.0008 

-0.0049 

2.9 

0.9999 

0.9969 

4.0782 

0.0146 

0.0004 

-0.0025 

u  0.01) 

0 

0 

0 

-0.0049 

1.3238 

10.0 

-7.8608 

0.2 

0.1559 

0.1323 

0.0218 

1.3227 

8.4406 

-7.8506 

0.4 

0.3108 

0.2640 

0.1016 

1.3076 

6.8922 

-7.7^ 

0.6 

0.4606 

0.3928 

0.2339 

1.2656 

5.3938 

-7.3084 

0.8 

0.5995 

0.5158 

0.4166 

1.1869 

4.0047 

-6.5846 

1.0 

0.7209 

0.8194 

0.6288 

0.6466 

1.1068 

2.7906 

-5.5467 

1.2 

0.7282 

0.9191 

0.9138 

1.6054 

-4.3042 

1.4 

0.8926 

0.8108 

1.2278 

0.7377 

1.0735 

-3.0373 

1.6 

0.9419 

0.8756 

1.5660 

0,5561 

0.5812 

-1.9262 

1.8 

0.9716 

0.9230 

1,9264 

0.3937 

0.2842 

-1.0920 

2.0 

0.9675 

0.9553 

2.3025 

0.2560 

0.1246 

-0.5464 

2.2 

0.9951 

0.9757 

2.6691 

0.1566 

0.0490 

-0.2431 

2.4 

0.9963 

0.9677 

3.0620 

0.0660 

0.0171 

-0.0949 

2.6 

0.9995 

0.99>^ 

3.4786 

0.0456 

0.0054 

-0.0326 

2.8 

0.9996 

0.9975 

3.8770 

0.0219 

0.0015 

-0.0096 

2.9 

0.9999 

0.9996 

4.0766 

0.0147 

0.0008 

-0.0051 

Table  4  (Contd.) 


T1 

u 

^0 

^2 

^3 

^4 

Bi 

(Ci^- 

0.2) 

0 

0 

0 

-0.1167 

1.2226 

0.2000 

-0.1521 

0.2 

0.1288 

0.1224 

-0.0876 

1.2488 

0.1742 

-0.1566 

0,4 

0.2644 

0,2471 

-0.0016 

1.2610 

0.1471 

-0.1588 

0.6 

0.4036 

0.3718 

0.1387 

1.2451 

0.1193 

-0.1559 

0.8 

0.5408 

0.4932 

0.3296 

1.1894 

0.0918 

-0.1456 

1.0 

0.66B3 

0.6070 

0.5664 

1.0886 

0.0663 

-0.1272 

1.2 

0.7778 

0.7087 

0.8435 

0.9464 

0.0444 

-0,1023 

1.4 

0.8636 

0.7948 

1.1545 

0.7757 

0.0273 

-0.0747 

1.6 

0.9239 

0.8633 

1.4929 

0.5957 

0.0152 

-0.0489 

1.8 

0.9617 

0.9142 

1.6526 

0.4265 

0.0076 

-0.0285 

2.0 

0.9627 

0.9495 

2.2275 

0.2837 

0.0034 

-0,0147 

2.2 

0.9930 

0,9721 

2.6130 

0.1749 

0.0014 

-0.0067 

2.4 

0.997*» 

0.9656 

3.0051 

0.0997 

0.0005 

-0.0027 

2.6 

0,9991 

0.9930 

3.i(011 

0.0526 

0.0002 

-0.0009 

2.8 

0.9996 

0.9^ 

3.7992 

0.0256 

-0.0003 

a.9 

0.9996 

0.9976 

3.9967 

0.0173 

0.0000 

-0.0002 

•  P-*) 


0 

0 

0 

-0.2968 

1.0561 

0.4000 

-0.2922 

0.2 

0.1000 

0.1071 

•0.2671 

1.1198 

0.3600 

-0.3110 

0.4 

0.2120 

0.2199 

-0.1742 

1.1706 

0.3152 

-0.3274 

0.6 

0.335*» 

0.3368 

-0.0246 

1.1948 

0.26^ 

-0.3357 

0.8 

0.4670 

0.4544 

0.1761 

1.1780 

0.2132 

-0.3295 

1.0 

0.5996 

0.5682 

0.4212 

1.1m 

0.1601 

-0.3039 

1.2 

0.7242 

0.6731 

0.7034 

0.9942 

0.1103 

-0.2587 

1.4 

0.8302 

0.7646 

1.0158 

0.8382 

0,0679 

-0,1999 

1.6 

0.9107 

0.8394 

1.3523 

0.6621 

0.0357 

-0,1363 

1.8 

0.9648 

0.89^ 

1.7077 

0.4879 

0.0141 

-0.0849 

2.0 

0.9967 

0.9376 

2.0774 

0,3343 

0.0013 

-0.0460 

2.2 

2.4 

1.0131 

1.0206 

0.9647 

0.9613 

2.4574 

2.8445 

0.2124 

0.1250 

-0.0052 

-0.0082 

-0.0219 

-0.0092 

2.6 

1.0235 

0.9907 

3.2360 

0.0680 

-0.0094 

-0.0034 

2.8 

1.0246 

0.9957 

3.6301 

0.0342 

-0.0098 

-0.0011 

2.9 

1.0248 

0.9971 

1 

_ i 

3.6^ 

0.0236 

-0.0099 

— 1 

-0.0006 
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TMble  U  (C«atd.) 


Ti 

WSmH 

u 

yo 

^2 

— 

^3 

^4 

(Cl,- 

0.6) 

0 

0 

0 

-0.6018 

— 

0.8110 

0.6000 

-0.3924 

0.2 

0.0666 

0.0846 

-0.5680 

0.9108 

0.5600 

-0.4306 

0.4 

0.1455 

0.1790 

-0.4685 

1.0095 

0.5127 

-0.4711 

0.6 

0.2386 

0.2823 

-0.3076 

1.09^ 

0.4568 

-0.5067 

0.8 

0.3464 

0.3919 

-0.0919 

1.1350 

0.3921 

-0.5279 

1.0 

0.4669 

0.5036 

0.1696 

1.1271 

0.3196 

1.2 

0.5937 

0.6120 

0.4670 

1.0590 

0.2438 

1.4 

0.7160 

0.7110 

0.7912 

0.935^ 

0.1704 

-0.4096 

1.6 

0.8216 

0.7960 

1.1350 

0.7728 

0.1070 

-0.3109 

1.8 

0.9017 

0.8642 

1.4935 

0.5950 

0.0590 

-0.2089 

2.0 

0.951^3 

0.9149 

1.8636 

0.4257 

0.0274 

-0.1233 

2.2 

0.9643 

0.9500 

2.2429 

0.2825 

0.0094 

-0.0638 

2.4 

0.9990 

0.9725 

2.6293 

0.1736 

0.0006 

-0.0289 

2.6 

1.0054 

0.9659 

3.0206 

0.0986 

-0.0032 

-0.0115 

2.8 

1.00^ 

0.9932 

3.4152 

O.05IB 

-0.0047 

-0.0040 

2.9 

1.0083 

0.9954 

3.6132 

0.0365 

-0.0050 

-0.0022 

<Clw 

-  0.8) 

0 

0 

0 

-1.2249 

0.3702 

0.8000 

-0.3993 

0.2 

0.0321 

0.0424 

-1.1983 

0.4733 

0.7743 

-0.4502 

0.4 

0.0713 

0.0960 

-1.1154 

0.7429 

-0.5096 

0.6 

0.1200 

0.1624 

-0.9730 

0.7368 

0.7040 

-0.57to 

0.8 

0.1807 

0.2424 

-0.7707 

0.8766 

0.6554 

-0.6451 

1.0 

0.2562 

0.3350 

-0.5127 

0.9957 

0.5950 

-0.7082 

1.2 

0.3486 

0.43^ 

-0.2076 

1.0695 

0.5211 

-0.7515 

1.4 

0.4578 

0.5429 

0.1321 

1.0779 

0.4337 

-0.7566 

1.6 

0.5790 

0.6463 

0.4936 

1.0139 

0.3368 

-O.707I 

1.8 

0.7013 

0.7404 

0.8662 

0.8867 

O.23B9 

-0.5990 

2.0 

0.8101 

0.8202 

1.2437 

0.7196 

0.1519 

-0.4502 

2.2 

0.8933 

0.8828 

1.6242 

0.5411 

0.0654 

-0.2959 

2.4 

0.9474 

0.9283 

2.0080 

0.3768 

0.0421 

-0.1690 

2.6 

0.9773 

0.9590 

2.3958 

0.2428 

0.0182 

-0.0840 

2.8 

0.9915 

0.9780 

2.7876 

0.1446 

0.0068 

-0.0364 

2.9 

0.9951 

_ 

0.9843 

2.9848 

0.1084 

0.0039 

-0.0228 
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